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' Abstract 

> 

. We study the motion of vortices in the planar Ginzburg-Landau model with 

Schrödinger-Chern-Simons dynamics. We compare the moduli space approximation 
with the results of numerical simulations of the füll field theory and find that there 
is agreement if the coupling constant is very close to the critical value separating 
Type I from Type II superconductors. However, there are significant qualitative dif- 
ferences even for modest deviations from the critically coupled regime. Radiation 
effects produce forces which are of the same order of magnitude as the intervortex 
force and therefore have a significant impact on vortex motion. We conclude that the 
moduli space approximation does not provide a good description of the dynamics in 
this regime. 
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1 Introduction 



The planar Ginzburg-Landau model provides a good mathematical description of the static 
properties of vortices in thin superconductors, but in order to describe vortex motion this 
needs to be supplemented by some appropriate dynamics. There are several possibilities 
for the form of this dynamics and it is currently difficult to determine which of these is the 
closest to modelling real superconductors, since in experimental situations the vortices are 
often pinned, so it is difficult to identify the underlying dynamical behaviour. 

It has been argued [1] that at very low temperatures dissipation can be ignored and 
vortex motion is orthogonal to the force acting, so that two vortices circulate around each 
other, as in a fluid. An interesting Lagrangian field theory to describe this type of behaviour 
has been proposed by Manton [4] in which the potential part is the usual Ginzburg-Landau 
energy and the kinetic terms are linear in the time derivatives of the complex scalar field 
and contain a Chern-Simons term. The equation of motion for the complex scalar field is 
a gauged Schrödinger equation, hence the name Schrödinger-Chern-Simons dynamics. It 
is Manton's model that we shall study in this paper. 

At critical coupling, which separates Type I from Type II superconductors, minimizers 
of the Ginzburg-Landau potential energy are static Solutions of the Schrödinger-Chern- 
Simons equations, so there is no vortex motion. Thus to investigate vortex dynamics the 
model needs to be studied away from critical coupling. 

The moduli space approximation [3] to soliton dynamics, in which the field theory 
dynamics is approximated by motion on a finite dimensional manifold, has been applied to 
Schrödinger-Chern-Simons vortex dynamics [4, 8]. The moduli space used is the manifold 
of static vortices at critical coupling, so it should be most accurate when the coupling 
constant is close to its critical value, though it can not equal this value if there is to be 
vortex motion, as mentioned above. In particular, it predicts that indeed two vortices 
circulate around each other at constant speed, and provides an expression for the period of 
this motion in terms of quantities defined on the moduli space. Polygonal configurations 
of iV vortices can also be studied and the moduli space approximation predicts that for 
well separated vortices these Symmetrie configurations are stable if and only if N < 6. 

In this paper we perform numerical simulations of the füll field theory dynamics and 
compare the results with the predictions of the moduli space approximation. For two 
vortices we find that there is a good agreement when the coupling constant is very close 
to the critical value. However, there are significant qualitative differences even for modest 
deviations away from critical coupling. We find that two vortices move on spirals, rather 
than circles, and spiral in or out depending upon the value of the coupling constant. This 
behaviour is due to radiation effects which are neglected in the moduli space approximation 
but turn out to produce forces which are of the same order of magnitude as the intervortex 
force. Polygonal arrangements also exhibit the same spiral phenomenon and the stability 
properties do not always agree with the moduli space predictions. 
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2 The Schrödinger-Chern-Simons Model 




The Ginzburg-Landau potential energy is given by 

ß 2 +Ä0A0+^(l-00) 2 ^ 2 X, (2.1) 

where is a complex scalar field, representing the electron pair condensate, and the co- 
variant derivative is = <9j0 — ia^, formed using the abelian gauge potential Oj with an 
associated magnetic field 

= ^02-0^01. (2.2) 

In (2.1), and the following, the summation Convention applies to the two spatial indices 
only, ie. i — 1,2. 

The finite energy topological solitons of this model are known as vortices. The vortex 
number N is a topological degree and is equal to the total magnetic flux in units of 2i\, 

N = ^- f B d 2 x. (2.3) 

2lt J 

N is also equal to the total number of zeros of the complex scalar field <fi counted with 
multiplicity. These zeros can be interpreted as the vortex positions when the vortices are 
well separated. 

The positive parameter A plays a crucial role in the properties of the model and its 
vortex Solutions. For A < 1 the superconductor is of Type I and the potential energy V of 
two vortices is an increasing function of their Separation ie. vortices attract. 

For A > 1 the superconductor is of Type II and the potential energy of two vortices is 
a decreasing function of their Separation ie. vortices repel. At the critical coupling A = 1, 
separating Type I from Type II superconductors, the potential energy is independent of 
their Separation. Moreover, at critical coupling, and for all positive vortex numbers N, the 
parameter space of all minimal energy N- vortices forms a 2iV-dimensional smooth manifold, 
known as the moduli space M^. In this case the potential energy of any N- vortex Solution 
is equal to rrN for arbitrary values of the vortex positions. 

At critical coupling the second order field equations which follow from the Variation of 
the potential (2.1) can be reduced to two first order Bogomolny equations 

D l( f) + iD 2 (f) = 0, (2.4) 
£>_!(!_ 00) = o, (2.5) 

and these will play an important role later. 

The Schrödinger-Chern-Simons model introduced by Manton [4] is a non-dissipative 
model for vortex motion in thin superconductors. The model is defined by the Lagrangian 

L = j Q (0A)0 - W) + (Ba + eia 2 - e 2 ai) - a ) d 2 x - V, (2.6) 
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where we have scaled some possible parameters to convenient values and is the electric 
field ti = dodi — di<i . Note that the term Ba + e\a 2 — e 2 ai is simply the Chern-Simons 
term written out explicitly and that the — üq term is required to allow the possibility of 
a condensate at infinity. The potential term V of the Lagrangian (2.6) is the Ginzburg- 
Landau energy (2.1). This Lagrangian is both gauge invariant and Galilean invariant, 
when one takes into account a constant external transport current. Although this is an 
important feature of the model we shall not be concerned with this aspect here. 
The equations of motions which follow from (2.6) are given by 

iD 0( ß = ^(1-00)0, (2.7) 

-EijdjB = % - (0A0 - 0Ä0) + 2e ij e j , (2.8) 

B = ^(1-00). (2.9) 

The first equation is a gauge covariant nonlinear Schrödinger equation. The second equa- 
tion is an Ampere equation, where the total current is the sum of the usual supercurrent 
and a Hall current orthogonal to the electric field. The third equation is a generalized 
Gauss law, which contains no time derivatives. It can be interpreted as a constraint on 
the initial data since it can be shown that the first two equations already imply that it is 
conserved, namely, 

|(fl-I(l-fc))=0. (2-10) 

Therefore, once equation (2.9) is satisfied for the initial data, equations (2.7) and (2.8) 
guarantee that it is satisfied for all later times. 

As the Lagrangian is linear in time derivatives then the kinetic energy makes no contri- 
bution to the conserved energy, which is simply the potential energy V, whose conservation 
is easily checked using the equations of motion. 

If the Schrödinger equation (2.7) and the Ampere equation (2.8) are considered for 
static fields, with a vanishing time component for the gauge potential üq = 0, then these 
two equations reduce to the two second order field equations obtained from the Variation 
of the Ginzburg-Landau energy (2.1). However, away from critical coupling, these static 
Ginzburg-Landau vortices will not satisfy the Gauss law (2.9) and hence are not static 
Solutions of the Schrödinger-Chern-Simons equations. The exception is at critical coupling 
A = 1 where Ginzburg-Landau vortices also satisfy the Bogomolny equations (2.4) and 
(2.5), the second of which is precisely the Gauss law (2.9). 

The first issue to consider is therefore the static Solutions of the Schrödinger-Chern- 
Simons model away from critical coupling. We shall restrict to axially Symmetrie Solutions, 
since it is expected that all static Solutions will have axial symmetry when A ^ 1. We work 
in the radial gauge a r = 0, with ag(r) and ao(r) funetions of the radius r only, and the 
complex scalar field having the Standard iV-vortex form = f(r)e tNe , where f(r) is the 
radial profile funetion. 
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Substituting this axial ansatz into the field equations (2.7), (2.8) and (2.9) yields the 
following set of ordinary differential equations 



f" + lf'-^(N-a e ) 2 f + ±(l-f)f + 2a f = 0, 
-a' e + f - 1 = 0, 



(2.13) 



(2.11) 



(2.12) 



where / denotes differentiation with respect to r. The boundary conditions for the fields 
are /(0) = 0, f(oo) = 1, a e (0) = 0, a e (oo) = N and a' (0) = 0. Note that from (2.11) these 
boundary conditions automatically imply that ao(oo) = 0. 

These equations are solved using a gradient flow method with a fictitious energy con- 
structed from the Square of the field equations. The results for the N — 1 vortex are 
displayed in Fig. 1, for three values of the coupling constant; A = 1 (solid curves), A = 0.5 
(dashed curves), A = 2 (dotted curves). The critical coupling case A = 1 is shown for 
comparison, where a = and the functions f(r) and a ö (r) are the usual fields of the 
Ginzburg-Landau vortex. The qualitative features of both the profile function f(r) and 
the angular component of the gauge field ag(r) do not vary significantly away from critical 
coupling. If A < 1 then /(r) is slightly wider and a g (r) is a little narrower, with the oppo- 
site being true for A > 1, that is, f(r) is narrower and ag(r) is wider. The most significant 
feature is the behaviour of the temporal component of the gauge potential a , which is 
positive for A < 1 and negative for A > 1. Thus, away from critical coupling, vortices not 
only have a magnetic field but also a tiny electric field. For a Single vortex the electic field 
is radial and positive if A < 1, whereas it is negative if A > 1. 

As we have mentioned earlier, the conserved energy in the Schrödinger-Chern-Simons 
model is the Ginzburg-Landau energy V. However, for A ^ 1 the static fields are not those 
of the Ginzburg-Landau model, so the energy dependence on A, for varying vortex numbers 
N, can not be inferred from knowledge of the Ginzburg-Landau model. In Fig. 2 we plot 
the energy per vortex (in units of ir) of the axial A^-vortex as a function of A for N — 1 
(solid curve), N = 2 (dashed curve), A^ = 3 (dotted curve). We have also computed similar 
results for larger values of N, but for clarity they are not shown on this plot. 

Despite the above comment we find that the qualitative behaviour of the energy is 
similar to Ginzburg-Landau vortices. In particular, for multi- vortices we see that for A < 1 
the energy of A" well separated vortices is higher than the energy of an axially Symmetrie 
charge N vortex, so vortices are attractive in this sense. For A > 1 vortices are repulsive, in 
that the energy of A" well separated vortices is lower than the charge A" axial configuration. 
At critical coupling, A = 1, the static vortices coincide with Ginzburg-Landau vortices, so 
the energy per vortex equals n for all N, as it does for all vortices in the 2A r -dimensional 
moduli space Mjv- 

The above features of the vortex interaction energy will play an important role in the 
dynamics of vortices, as we shall see later. 
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Figure 1: Graphs of the functions f(r), a 9 (r) and a (r) for couplings A = 1 (solid curves), 
A = 0.5 (dashed curves), A = 2 (dotted curves). 



3 The Moduli Space Approximation 

The main idea of the moduli space approximation [3] is to project the dynamics of the füll 
field theory onto a suitable finite dimensional space of field configurations. In the simplest 
Situation, such as vortices at critical coupling or BPS monopoles, this finite dimensional 
space is the moduli space of static minimal energy Solutions, with a given soliton number. 
This approach is well established (for a review see [6]) for relativistic Lagrangians, where 
the dynamics is second order in time. Restricting the field theory Lagrangian to the moduli 
space produces (upto an irrelevant constant) a purely kinetic Lagrangian on the moduli 
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Figure 2: The energy per vortex (in units of n) as a function of the coupling A for the axial 
iV-vortex with N — 1 (solid curve), N = 2 (dased curve), N = 3 (dotted curve). 



space, which may be interpreted as the moduli space metric. The slow motion of solitons 
is approximated by geodesic motion with respect to this metric. For relativistic dynamics 
of Ginzburg-Landau vortices in the abelian Higgs model the validity of this approximation 
has been proved with mathematical rigour [13]. 

A slightly more complicated Situation arises if a theory is considered at parameter values 
for which the moduli space of static Solutions is not large enough to describe solitons with 
arbitrary positions, for example, vortices in the abelian Higgs model away from critical 
coupling. In such situations the dynamics can be truncated to motion on the moduli space 
of the theory at critical coupling, but now there will also be a potential energy function 
on the moduli space, so the dynamics is no longer described by geodesic motion. This 
approach has been applied in different models [14, 10], and there is a good agreement with 
results obtained from numerical simulations of the füll relativistic field theory dynamics. 
For relativistic vortex dynamics near to critical coupling rigourous error estimates can 
again be obtained [13]. 

Schrödinger-Chern-Simons vortex dynamics is first order in time, so initial conditions 
consist only of vortex positions. Furthermore, at critical coupling the Ginzburg-Landau 
vortices are static Solutions, so there is no vortex motion. Thus, in constrast to relativistic 
dynamics, it is necessary to move away from critical coupling in order for there to be 
any vortex dynamics to study. Manton has argued [4] that it is possible to use moduli 
space techniques to describe Schrödinger-Chern-Simons vortex dynamics close to critical 
coupling, and this has been further investigated in some detail by Romäo and Speight[8]. 
We summarise the results of the relevant calculations below. 

Let q a , for a = 1, . . . , 2N, be real coordinates on the moduli space M N) for example the 
real and imaginary parts of the vortex positions when space is identified with the complex 
plane. As the Schrödinger-Chern-Simons Lagrangian (2.6) depends only linearly on time 
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derivatives then the Lagrangian restricted to the moduli space must have the form 

L Ted = A a (q)q a - V red (q), (3.1) 

where A a (q) may be interpreted as the components of a U(l) connection form A = A a dq a 
on Mtv- The first order equations of motion which follow from (3.1) are 

= (3 ' 2) 

where T = dA is the curvature form of A. In order to study the details of this reduced 
dynamics we first need a better understanding of the moduli space M N . 

By considering the Bogomolny equations (2.4) and (2.5) Taubes [15] showed that it is 
possible to eliminate the gauge potential as follows. Define 

h = log\<j>\ 2 , (3.3) 

which is gauge invariant and finite, except at the zeros of <fi. Using the first Bogomolny 
equation (2.4), it can be shown that the magnetic field can be expressed as 

B = ~\^ 2 h. (3.4) 

The second Bogomolny equation then becomes 

V 2 h + 1 - e h = 0. (3.5) 

Recall h has logarithmic singularities at the zeros of <fi. This can be accounted for by 
including delta-function sources into the equation, 

N 

V 2 h + 1 - e h = An £ 6(z - Z r ), (3.6) 

r=l 

where 5 is the two dimensional delta-function and Z r , with r — 1, . . . , N., are the comp lex 
positions of the zeros of <fi in R 2 , which is identified with the complex plane. It can be 
shown that h has the following expansion around a zero z = Z r 

h{z,z) = 2 log \z - Z r \ + a r + ]-b r (z - Z r ) + h r (z - Z r ) + 0(\z - Z r \ 2 ). (3.7) 

It turns out that the coeflicients b r contain all the necessary Information to reconstruct the 
geometry of the moduli space. 

Following the work of Strachan [12] , Samols [9] was able to show that the moduli space 
metric derived from vortex dynamics in the relativistic abelian Higgs model is given by 



N ( db s \ - 

ds 2 = 7r \örs + 2 -q^-J dZ r dZ s , (3. 
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which is a Kahler metric. The corresponding Kahler form is 

N / QI \ 

w = yE [8rs + Igf) dZ r A dZ s . (3.9) 

and it turns out that the curvature form J 7 , which appears in Schrödinger-Chern-Simons 
vortex dynamics (3.2), is simply proportional to this Kahler form, namely T = —2uj. 

Note that the above formulae for the geometrical objects on the moduli space are rather 
implicit, since they involve the expansion coemcients b r , and it is not possible to determine 
these explicitly in terms of the vortex positions Z r . Furthermore, the reduced potential 
^red' wn i cn i s the other ingredient required in the vortex equations of motion (3.2), is also 
not available explicitly. It is given by 

V ied = Nn + ^J(l- |0| 2 ) 2 d 2 *, (3.10) 

where the field in the above equation is the Solution of the Bogomolny equations with 
zeros at Z r . 

Despite the fact that the quantities which appear in the moduli space approximation are 
not known explicitly it is still possible to derive some interesting results. As an example, 
it is possible to show that the naive centre of the iV-vortex system 

1 N 

Z=~Y,Zr (3.11) 

JV r=l 

is a constant of the motion. 

The qualitative motion of two vortices can also be understood. As we have just 
mentioned, the centre Z is a constant of motion, which we may choose to be the ori- 
gin. The vortex positions can then be written as Z\ = —Z 2 = pe' ld , and it follows that 
bi = —62 = b(p)e td . Furthermore, the reduced potential V^. e( j depends only on p. In this 
case the equations of motion (3.2) become 

p = 0, (3.12) 



2^ (p 2 + 2p b{p)) 9 = --^Spl. (3.13) 

These equations imply that both the vortex Separation 2p and the angular velocity 9 are 
constant. So, the vortices circle around each other at constant speed and Separation. If 
A < 1 then dV ie ^/dp is positive, so the two vortices circle each other clockwise, whereas 
for A > 1 the motion is anticlockwise. 

An expression for the period T of this circular motion is given by 

^l^ + 2pb(p)) 



T = 



dV ie ^/dp 

Note that the combination |1 — A|T is a function only of p, that is, it is independent of A, 
so the A dependence of the period is known. If p is large then the speed is exponentially 
small. For small p then both the numerator and the denominator in (3.14) are 0(p 3 ), so 
the period approaches a finite non-zero limit as p — > 0. 
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4 Two Vortex Dynamics 



In the remainder of this paper we compare the detailed predictions of the moduli space 
approximation with the results of numerical simulations of the füll field theory dynamics. 
We begin in this section with a study of the motion of two vortices. 

There are no explicit expressions for the functions b(p) and V rec {(p) which appear in 
the 2-vortex period formula (3.14), but they can be computed numerically following the 
approaches described in [9] and [10]. This involves numerically solving equation (3.6) to 
compute h, for which we use a gradient flow algorithm with a change of variables that 
move the delta function source terms to infmity. The Solution h is calculated for a ränge 
of values for the vortex Separation s = 2p, and this is used both to read-off the function 
b(p) using the definition (3.7), and to calculate the reduced potential (3.10) by performing 
a numerical Integration over the plane, using the simple relation (3.3) between h and |0| 2 . 

As mentioned in the previous section, the combination |1 — X\T, which we refer to as 
the scaled period, is independent of A. It is this quantity that is plotted in Fig. 3, as a 
function of the Separation s = 2p, for 1.5 < s < 3.6. This graph displays the expected 
feature that the period of rotation increases with the vortex Separation. 

Recall that the period tends to a finite non-zero value as the Separation tends to zero, 
and it would be useful if our numerical computation could be extended to the ränge < 
s < 1.5, so that this value could be computed. However, recall that the period is calculated 
as the ratio of two quantities which are both 0(p 3 ) as p tends to zero. Furthermore, both 
of these quantities involve differentiating terms which we compute numerically. Thus, 
although the numerically computed functions b(p) and V re( i(p) are reasonably accurate, 
calculating the period for small separations is much more sensitive to numerical errors, 
since it involves the ratio of two terms whose quartic approaches to zero need to be found. 
This is why our calculation for the period is only presented for p > 0.75, since for smaller 
values we are not able to achieve the accuracy required. To improve the accuracy the 
Solution h of the partial differential equation (3.6) would need to be solved on larger grids 
than we currently use, but this is computationally rather expensive. 

An alternative approach would be to linearize around the axially Symmetrie two vortex 
Solution in order to obtain the period of rotation for minimal Separation. Since the two 
vortex Solution is not known analytically, this is a nontrivial numerical problem, but it is 
an interesting project for further study. 

We now turn to numerical simulations of the füll field theory dynamics. The evolution 
equations (2.7) and (2.8) are solved using a second order finite difference scheme for the 
spatial derivatives with lattice spacing dx = 0.1 on a 400 x 400 grid. The time evolution 
uses a fourth order Runge-Kutta algorithm with a time step dt = 0.005. Note that the 
time step is much smaller than the lattice spacing, as expected for the numerical Solution 
of a partial differential equation which is first order in time derivatives and second order 
in spatial derivatives. 

Recall that the Gauss law equation (2.9) is a constraint on the initial data. This means 
that it is a non-trivial task to construet initial conditions for simulations; in particular the 
usual method of a simple produet ansatz for well separated vortices can not be used since 
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Figure 3: The scaled period T|l — A| as a function of the vortex Separation s, as predicted 
by the moduli space approximation. The data points show the scaled period as calculated 
from füll field theory simulations with couplings A = 0.9 (circles), A = 0.75 (squares), 
A = 0.5 (triangles). 
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Figure 4: The vortex positions for two vortices with initial Separation s = 2.5 and coupling 
A = 0.9. 
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this will not satisfy the constraint. One of the main motivations for our work is to compare 
field theory dynamics with the moduli space approximation, and this suggests a natural 
Solution to this problem. As initial conditions we use Solutions of the Bogomolny equations 
(2.4) and (2.5), which are static Solutions only at critical coupling, since these automatically 
satisfy the Gauss law constraint (2.9). In summary, to obtain initial conditions for our 
simulations we fix the vortex positions, then solve equation (3.6) to obtain the field h with 
these vortex positions, and finally reconstruct the initial condition fields <fi and from h, 
setting the initial value of a to zero. 

The results of a two vortex Simulation with A = 0.9 and initial positions (xi,x 2 ) = 
(±1.25,0), so that the initial Separation is s = 2.5, are presented in Fig. 4. The vortex 
positions are defined as the locations of the zeros of the Higgs field, and it is these that 
are shown in Fig. 4 for times < t < 1646. The two vortices circle around each other 
clockwise at approximately constant speed and Separation. The results are displayed for 
a half-period giving T/2 = 1646, which corresponds to more than 300,000 timesteps with 
dt = 0.005. The motion is very slow because the coupling A = 0.9 is near to the critical 
value, so the forces are weak. 

Clearly this Simulation confirms the qualitative prediction of the moduli space approxi- 
mation, that two vortices circle around each other. To compare the quantitative prediction 
we plot the period of this motion against the Separation as a circle data point on Fig. 3. 
The Separation is slightly larger than 2.5 since we calculate it not from the initial value, 
but as a value averaged over the half-period, and the vortices very slightly drift apart. The 
additional circles in Fig. 3 represent the results of similar computations to calculate the 
period for different values of the initial Separation. The circles in Fig. 3 are reasonably 
close to the solid curve prediction of the moduli space approximation, indicating that at 
this value of the coupling it provides good results. 

In Fig. 5 we present the results of a Simulation with initial Separation s = 2 and coupling 
A = 1.5. For clarity we display the position of only one of the vortices; the position of 
the other vortex being obtained by symmetry as {x\,x 2 ) i— > (— x±,— x 2 ). Fig. 5(a) shows 
the position in the (xi,x 2 ) plane for times < t < 3670, and Fig. 5(b) displays the X\ 
component of the position as a function of time. 

From these figures it can be seen that this time the vortices circle around each other 
anticlockwise, as predicted by the moduli space approximation when A > 1. However, the 
Separation is far from being constant and increases substantially with time, producing an 
outward spiral motion of the vortices. The outward motion of the vortices is the same order 
of magnitude as the circular motion, indicating that there are significant forces acting on 
the vortices that are neglected in the moduli space approximation. Fig. 5(b) displays the 
X\ component of the position as a function of time, verifying that the time taken for the 
vortices to circle each other increases as they move further apart. 

A natural Interpretation of the outward spiral motion is that the vortices radiate as 
they circulate. Recall from Section 2 that for A > 1 the energy of two vortices decreases 
with increasing Separation, so an energy loss due to radiation would result in the vortices 
moving away from each other. If this is indeed the correct explanation of the outward 
spiral motion, then it predicts that for couplings not too close to the critical value and 
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Figure 5: The position of one of the vortices, for times < t < 3670, for a two vortex 
Simulation with initial Separation s = 2 and coupling A = 1.5. (a) The position in the 
(xi,X2) plane, (b) The x\ component of the position as a function of time. 



satisfying A < 1 then the vortices should circle around each other with an inward spiral 
motion. This is because for A < 1 the results from Section 2 show that the energy of two 
vortices decreases as the Separation decreases, so this time an energy loss due to radiation 
would result in the vortices moving towards each other. In order to test this prediction we 
perform a Simulation with initial Separation s = 7 and coupling A = 0.5. The results are 
displayed in Fig. 6 using the same format as in Fig. 5. 

The vortices circle each other clockwise and the inward spiral motion is clearly visible 
in Fig. 6. This is consistent with the above expectation based on an energy loss due to 
radiation. The Xi component of the position, plotted in Fig. 6(b), suggests that the energy 
loss will eventually lead to the final configuration being a static N = 2 axial vortex, though 
the approach to this static Solution is very slow. The Simulation presented already involves 
one million time steps of the evolution algorithm, so it is computationally too expensive 
to attempt to follow the evolution for any substantially longer time. There is always the 
possibility that the configuration could stabilize to a rotation at a very small non-zero 
Separation, but this seems unlikely. 

It is evident from Fig. 6(a) that for later times there is a significant wobbling of the 
vortex motion. This wobbling component eventually appears in all our simulations. It 
can be seen in Fig. 4 and Fig. 5, though it is much less pronounced in these cases. We 
do not know whether this is a purely numerical artifact or has some physical significance. 
It certainly signals the onset of a numerical instability in our evolution algorithm since 
our codes are unable to continue the evolution for very long after a pronounced wobbling 
emerges. Monitoring the constraint (2.9) we find that it is satisfied to within a reason- 
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Figure 6: The position of one of the vortices, for times < t < 5000, for a two vortex 
Simulation with initial Separation s = 7 and coupling A = 0.5. (a) The position in the 
(xi,X2) plane, (b) The x\ component of the position as a function of time. 



able numerical accuracy until shortly before the evolution algorithm fails, but we do not 
know whether a drift away from this constraint is a cause or Symptom of the instability. 
There is an initial release of radiation at the start of all our simulations, since the initial 
conditions are created for vortices with critical coupling whereas the simulations require a 
non-critical value. There is some evidence, from changing grid sizes, that the interaction 
of this radiation with the boundary of the grid may be a factor in limiting our Simulation 
times. We use Neumann boundary conditions, which results in some radiation reflecting 
from the boundary. It would be better if absorbing boundary conditions could be used, but 
this is not easy to implement since radiation in the Schrödinger-Chern-Simons equations 
is non-standard. Linearizing the evolution equations shows that the plane wave Solutions 
are elliptically polarized. See the appendix for further details. 

In summary, we have seen that for two vortices the moduli space prediction of circular 
motion of the vortices around each other is only a good approximation very close to critical 
coupling. For modest deviations from the critically coupled regime there is a substantial 
spiral motion, which is an inward spiral for Type I vortices and an outward spiral for 
Type II. Very close to critical coupling we found that the moduli space approximation for 
the period of the circular motion was quite accurate. Further away from critical coupling 
the motion is not even approximately periodic, but an indication of the time scale of the 
rotation can be computed by calculating the 'period' as twice the time taken for the vortices 
to make a rotation through 180°. We include values of these scaled periods in Fig. 3 for 
various separations and couplings A = 0.5 (triangles) and A = 0.75 (squares). Comparing 
these values with the near critical value A = 0.9 (circles) suggests that the moduli space 
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prediction for the period is indeed approached as the coupling tends to the critical value 
A = 1. 

5 Rotating Polygons 

Symmetry of the moduli space dynamics (3.2) implies that if there are N vortices at the 
vertices of a regulär iV-gon then this polygon of vortices will rigidly rotate about its centre 
with constant size and speed [4]. 

The linear stability of such a rotating polygon has been studied within the moduli space 
approximation [8]. This Situation is a gener alization of the two vortex system, so it is not 
surprising that neither the Kahler form nor the reduced potential are explicitly known. 
The approach of Ref. [8] is to consider the limit in which the size of the polygon is large, 
so that the vortices are well separated. For well separated vortices there is an asymptotic 
formula for the Kahler form [5] and the reduced potential can be approximated by the 
Ginzburg-Landau energy, for which there is an asymptotic formula in terms of the sum 
over the asymptotic two- vortex interaction energies [11]. With all these approximations, 
the result [8] is that the iV-gon is stable if and only if N < 6. 

Given the field theory simulations described in the previous Section, together with 
the energy calculations of Section 2, we expect that Symmetrie polygons will spiral in for 
A < 1 and out for A > 1, so in that sense they are all unstable. However, the stability 
results referred to above still apply if we ask about the stability of the symmetry of the 
configuration, that is, we shall refer to the Cn Symmetrie polygon as being stable if the 
Cn symmetry is stable, irrespective of a size instability. 

In this Section we present the results of füll field theory simulations on rotating poly- 
gons, to compare with the stability predictions of the moduli space approximation. These 
predictions are most likely to hold close to critical coupling and for polygons of a large size, 
so that the vortices are well separated. Unfortunately, both these regimes are difncult to 
simulate numerically since the forces are very weak. This makes the vortex motion very 
slow and requires an unacceptably long Simulation time. Nonetheless, we find that there 
is a substantial agreement between the moduli space predictions and the results of our 
simulations, though there is a disagreement for a pentagonal arrangement. 

The simulations we present are with a coupling constant A = 0.5, though we also 
obtained similar results with other values. There is no need to explicitly break the cyclic 
symmetry of the vortex positions, since the numerical discretization, and in particular the 
boundary of the grid, are sufneient to allow the symmetry to break if there is an unstable 
mode. 

In Fig. 7(a) we display the vortex positions for an initial triangulär arrangement with 
three vortices on a circle of radius 4. The white circles denote the initial positions and 
the lines indicate the resulting motion for times t < 2000. For clarity we do not display 
the vortex positions as lines for later times, but do show the positions (black circles) at 
the end of the Simulation at time t = 5000. At the end of the Simulation the vortices have 
made just over one füll rotation, and still preserve the triangulär symmetry. Note that at 
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Figure 7: Vortex positions for polygonal initial conditions. The white circles indicate the 
initial vortex positions and the black circles indicate the final vortex positions at the end 
of the Simulation. The lines show the vortex positions in the early stages of the Simulation, 
(a) A triangulär arrangement with vortices initially on a circle of radius 4. Lines are drawn 
for times t < 2000 and the black circles are the positions at the time t = 5000. (b) A 
pentagonal arrangement with vortices initially on a circle of radius 5. Lines are drawn for 
times t < 2500 and the black circles are the positions at the time t = 3200. 



t = 2000 the vortices have made less than one third of a füll rotation, but at t — 5000 
they have made just over a füll rotation. This is because the vortices radiate and their 
Separation decreases, with a corresponding decrease in the rotation period. This Simulation 
is consistent with the asymptotic moduli space prediction that a triangulär arrangement 
is stable. 

An arrangement of four vortices on the vertices of a Square produces a similar result, 
with the vortices moving on an inward spiral but preserving the C4 symmetry. 

The results of a Simulation with five vortices on the vertices of a regulär pentagon are 
shown in Fig. 7(b). The initial vortex positions (white circles) lie on a circle of radius 5, and 
the lines track the motion for t < 2500. It can already be seen from these vortex trajectories 




Figure 8: Energy density plots for an initial pentagonal arrangement of 5 vortices. The 
plots are at times t — and t = 2500. 
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Figure 9: Energy density plots for an initial hexagonal arrangement of 6 vortices. The 
plots are at times t — and t = 2000. 

that the cyclic C 5 symmetry gets broken, with one of the vortices already moving closer 
to the origin than the other four. The black circles denote the vortex positions at time 
t = 3200, and the pentagonal symmetry has obviously been destroyed, with one of the 
vortices being much further from the origin than the other four. The broken symmetry 
is most easily visualized through energy density plots. Fig. 8 shows energy density plots 
at times t — and t = 2500, clearly displaying the broken C 5 symmetry. Note that the 
vortices are more localized in the initial conditions than at later times. This is because 
the initial conditions are formed from critically coupled vortices, but these quickly relax 
to vortices with a localization appropriate for the chosen coupling A = 0.5. The numerical 
grid extends far beyond the inner portion plotted in these figures, but for clarity we do not 
display the entire grid. 

Further evolution of this configuration suggests that the vortices are tending towards 
the formation of a pair of axial N = 2 vortices, plus a Single N — 1 vortex, though our 
simulations can not be run long enough to confirm the final state. A reasonable expectation 
is that the vortices combine, at first in pairs, and eventually radiation loss will lead to a 
final configuration of an axial 5-vortex. 

The instability of the pentagonal arrangement disagrees with the asymptotic moduli 
space prediction and therefore raises a number of interesting questions as to the source of 
the disagreement. One possibility is that the asymptotic approximation is correct but that 
the Simulation is not in the asymptotic regime. This would suggest that there may be a 
critical size for the pentagonal arrangement, with a stability only if the pentagon is above 
the critical size. This could be investigated by numerically Computing the appropriate 
Kahler form and reduced potential, and then numerically solving an eigenvalue problem 
to determine if there are any negative modes. This is an interesting project for future 
investigation. 

Another possibility is that there is a stability region in coupling space, so that the 
pentagon is stable only if A is sufnciently close to critical coupling. It might be possible to 
investigate this within a moduli space approach by including the deformations of the vortex 
fields, such as the tiny electric field, which arise away from critical coupling and are not 
captured by restricting to the critically coupled moduli space. For A > 1 an appropriate 
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moduli space could be constructed from the unstable manifold of the axial N = 5 vortex, 
in a similar way that has been used [2] to study Schrödinger dynamics of two ungauged 
Ginzburg-Landau vortices. 

A final possibility is that the instability is due to a phenomenon that can not be 
captured by any Standard moduli space approach, and the interaction between vortices 
and radiation plays a crucial role. 

The asymptotic moduli space approach predicts that 6 vortices on the vertices of a 
regulär hexagon will be unstable. The results of such an initial arrangement, with the 
vortices on a circle of radius 5, are presented in Fig. 9. This figure displays energy density 
plots at times t — and t = 2000. Clearly the cyclic C 6 symmetry is broken by this time, 
so this result is in agreement with the instability prediction. Note from the second plot in 
Fig. 9 that the Cq symmetry is broken to a Ci symmetry at this stage in the evolution, 
and again there appears to be a tendency for vortices to form pairs. 

In summary, we find that the stability properties of regulär polygonal arrangements of 
vortices is in reasonable agreement with the moduli space prediction, in that we find that 
an iV-gon is stable only for small enough values N < N c . However, we conclude that the 
critical value is N c = 5 rather than N c = 6, and have suggested some approaches that 
could be used to investigate this discrepancy further. 



6 Hyperbolic Vortices 

The equations for static critically coupled Ginzburg-Landau vortices on the hyperbolic 
plane of curvature — \ are integrable, and this allowed Strachan [12] to calculate a general 
formula for the Kahler form on the moduli space M^. In this respect hyperbolic vortices are 
therefore simpler than Euclidean vortices. In this Section we exploit this fact to study the 
moduli space approximation to Schrödinger-Chern-Simons vortex dynamics on hyperbolic 
space. 

Consider a two-dimensional Riemannian manifold X, with metric 

ds 2 = Q(xi, X2)(dx\ + dx\). (6.1) 
Ginzburg-Landau vortices on X are minima of the energy 

V = \j x (tt^B 2 + DrfDrf + DrfD 2 (/> + ^(1 - 00) 2 ^) d 2 x. (6.2) 

At the critical coupling A = 1 the second order equations can again be reduced to first 
order Bogomolny equations, as in the Euclidean case. Introducing the function h, as in 
(3.3), then the gener alization of the Euclidean equation (3.6) becomes 

N 

V 2 h + tt-tte h = 4vr^5(a;-X r ), (6.3) 

r=l 

where V 2 is the Standard flat space Laplacian and X r G X are the vortex positions. 



18 



In the Poincare disc model, the metric of the hyperbolic plane with curvature — \ is 

g 

ds 2 = 7 . .„.„ dzdz = f2 dzdz, (6.4) 

(l — \z\ l Y 

where \z\ < 1. Setting h = 2g + 21og|(l — \z\ 2 ) the equation for h becomes Liouville's 
equation with sources 

TV 

V 2 g-e 2g = 27iY / S(z-Z r ), (6.5) 

r=l 

which can be solved exactly. The Solution is 

2 

(6.6) 



£ = -logi(l-m 2 ) + ^ 



df 



dz 



where f(z) is an arbitrary, complex analytic function. With a simple choice of phase the 
scalar field is given by 

i - Nf df 

* = T^TfTT; (67) 

Note that <p vanishes at the zeros of f. To ensure that has N zeros, is nonsingular inside 
the disc \z\ < 1 and is zero on the boundary \z\ = 1, then the function f(z) needs to have 
the Blaschke product form 



K+1 - z-a, 



fiz) = n fri 

1=1 Vi u 



(6.8) 



where |aj| < 1. 

Let us now specialize to the case of N vortices on an iV-gon, including the degenerate 
case N = 2. Such a configuration is described by the holomorphic function 

z (z N - a N ) 
1 - ä N z J 



The positions of the vortices are given by the zeros of the derivative df /dz, and these are 
the regulär iV-gon vertices given by 

z = a e 2mk/N (6.10) 
for k — 1, . . . , N. The parameter a can be expressed in terms of a via a = cry where 





a 


2Nf + 4N 2 


a 


2N _ l _ 


a 


2N 


2 


OL 


2N (N- 


■1) 







Using the results of [12] the Kahler metric on this two-dimensional submanifold of M/v can 
be calculated explicitly and is given by 

ds 2 = ^NM 2N ~ 2 d* da ( 2iV(l + |q| 2Af ) \ 

(1-|«| 2JV ) 2 V ^ (l - \a\™) 2 + AN*\a\™ ) ' 
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Note that for N = 2 this metric becomes the one derived by Strachan [12], which he used 
to study the relativistic dynamics of two hyperbolic vortices. 

Given the metric (6.12) the Kahler form u can simply be read off. Let a = re 10 , then 

u = f( r )dr Ad6, (6.13) 

where 




27riV 3 2N-i ( 2N(l + r 2N ) 



[l_ r 2N f \ _ r 2N^ + 4N 2 r 2N 

Now that we have the Kahler form, all that we require to study Schrödinger-Chern- 
Simons vortex dynamics is the reduced potential. This is given by the expression 

V ied = Nn + ^f x 0- ~ I0I 2 ) 2 « d 2 x, (6.15) 

where the explicit Solution for <fi is known from above using (6.7) and (6.9). It can be 
shown that this reduced potential depends only on \a\. Unfortunately this integral over 
the hyperbolic plane can not be performed explicitly (except for the special value a — 0) 
even in the case N = 2. The angular part of the integral can be calculated explicitly but 
the radial Integration must then be computed numerically. 

With the Kahler form and reduced potential given above the moduli space dynamics 
reduces to the equations of motion r = and 

2f( r )9 = (6.16) 

Thus, as in the Euclidean case, the configuration rigidly rotates with no change in size and 
with constant angular velocity 9. The period is 



T 



4vr/(r) 



dV Ted /dr 



(6.17) 



In order to compare the period calculation in hyperbolic space to the earlier results in 
Euclidean space, we need to calculate the proper distance. From the Poincare disc coor- 
dinate z = \z\e %x we make the transformation p = 2^/2t&nh~ 1 (\z\) so that the hyperbolic 
metric becomes 

ds 2 = dp 2 + 2 sinh 2 (j^J d X 2 - (6.18) 

Thus p should be equated with the radial distance in Euclidean space. 

In Fig. 10 we plot (dashed curve) the scaled period T\X — 1| as a function of the Separa- 
tion s = 2p for the two-vortex case ie. N = 2. For comparison we also reproduce the earlier 
Euclidean result (solid curve). The Euclidean and hyperbolic results are qualitatively sim- 
ilar. Both periods are monotonically increasing functions of the Separation and tend to 
finite non-zero limits as the Separation tends to zero. The simplifications of hyperbolic 
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Figure 10: The scaled period T\ \ — 1| against Separation s for two vortices in Euclidean 
space (solid curve) and hyperbolic space (dashed curve). 

space allow the period to be calculated accurately all the way down to zero Separation, 
and this value appears similar to the one extrapolated from the Euclidean data. Note that 
it is not expected that the Euclidean and hyperbolic periods agree at zero Separation since 
the vortices are extended objects, so even when both vortices are located at the origin they 
still feel the curvature of hyperbolic space. 

The period for an iV-gon arrangement of hyperbolic vortices can also be calculated 
using the above methods and the results are qualitatively similar to the 2-vortex case. The 
period at zero Separation is not very sensitive to the value of N but does slightly decrease 
as N increases. 

The similarity between the moduli space results for Euclidean and hyperbolic space 
suggests that it would be interesting to investigate the stability of polygonal arrangements 
of hyperbolic vortices. Although this is still not an easy exercise it is substantially simpler 
than the Euclidean case, since the exact static vortex Solutions are available at critical 
coupling. 

7 Conclusions 

In this paper we have studied the dynamics of vortices in Manton's Schrödinger-Chern- 
Simons model, and compared the results of füll field simulations with predictions from 
the moduli space approximation. We found that there is a good agreement for couplings 
extremely close to the critical coupling, but away from this value there are significant 
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qualitative differences, which we attribute to radiation effects which are not captured by 
the moduli space approach. This is of physical relevance since real superconductors are 
generally not close to the Type I/Type II transition. 

A novel Suggestion might be to try and modify the non-dissipative moduli space dy- 
namics to a new dissipative dynamics on the moduli space that would capture the energy 
loss effects. The Schrödinger-Chern-Simons flow is orthogonal, in a rigourous mathemat- 
ical sense, to the gradient flow of the Ginzburg-Landau energy, so a flow which contains 
both components should certainly reproduce the qualitative features that we have found. 
However, it is not clear how to derive such a flow from the field theory, but it would seem 
a worthwhile avenue for further investigation. 

In the Schrödinger dynamics of two ungauged Ginzburg-Landau vortices it has been 
shown analytically that two vortices radiate while rotating around each other and the 
asymptotic rate at which they radiate has been derived [7]. In the ungauged System the 
coupling constant plays no role. Motivated by the numerical results we have presented here, 
it would be interesting if the methods used in [7] could be applied to the Schrödinger-Chern- 
Simons model and similar results derived, including the coupling constant dependence. 
This is another interesting topic for future work, though the inclusion of gauge fields, and 
in particular the requirement of working with gauge independent quantities, appears to 
make the problem significantly more complicated than in the ungauged Situation. 

Finally, we have demonstrated that, at least within the moduli space approximation, 
the qualitative features of Euclidean vortices appear to be shared by hyperbolic vortices. 
This opens up the possibility of future studies on the stability properties of Symmetrie 
arrangements of hyperbolic vortices, for arbitrary separations, exploiting the fact that 
hyperbolic vortices are simpler to study due to the existence of exaet static Solutions at 
critical coupling. 



Appendix: The linearized field equations 

In this appendix, we discuss the linearized field equations and their plane wave Solutions. 
It is convenient to work in a gauge such that is real. Setting = 1 + er, we can linearize 
around the vaeuum = 1 and a ß = 0, by expanding the field equations (2.7), (2.8) and 
(2.9) in a and a ß to linear order. First, we eliminate a using the linearized real part of 
(2.7) 

2a = -VV + Xa, (A.l) 
where V 2 = df + d\. Then we replace a using the linearized Gauss equation 

a = -B. (A.2) 

Finally, from (2.8) we obtain 




(1 + A - V 2 )<9i<9 2 1-(1 + A- V 2 )<9 2 
-1 + (1 + A - V 2 )<9 2 -(1 + A - V 2 )9i9 2 




(A3) 
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Note that the determinant of the symbol of the differential Operator in (A.3) vanishes. In 
order to find plane wave Solutions we make the ansatz = Ajexp(i(k • x — out)). This 
results in a linear homogeneous matrix equation which has nontrival Solutions only if the 
determinant vanishes, namely, 

-Au 2 + k 4 + (1 + X)k 2 + 1 = 0, (A.4) 

where k 2 = k\ + k\. This gives rise to the nonlinear dispersion relation 

u = ^k 4 + (1 + A)fc 2 + 1, (A.5) 

which for critical coupling A = 1 reduces to 

uj = \{k 2 + l). (A.6) 

For general A the corresponding linearized Solution is given by 

cos (k • x — out + 5) \ 

ggg cos ( k • x ~ + S) + sin (k ■ x - ut + 5) ) ' 

(A.7) 

where et and 5 are constants. This shows that the plane wave Solutions are elliptically 
polarised, as can be seen for example from the fact that the modulus of the vector (A.7) 
is not time independent. 
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